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The most straightforward manner of determining masses and radii of neutron 
stars is by measuring the gravitational redshift of spectral lines produced in the 
neutron star photosphere; such a measurement would provide direct constraints 
on the mass-to-radius ratio of the neutron star, and therefore on the equation of 
state for neutron star matter. Using data taken with the Reflection Grating 
Spectrometer on board the XMlVI-yVewtow observatory we identify, for the first 
time, significant absorption lines in the spectra of 28 bursts of the low-mass X-ray 
binary EXO 0748-676. The most significant features are consistent with the Fe 
XXVI and XXV n = 2-3 and O VIII n = 1-2 transitions, with a redshift of z = 0.35, 
identical within small uncertainties for the diRerent transitions. This constitutes 
the first direct and unambiguous measurement of the gravitational redshift in a 


neutron star. 
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The XMM-A^^Hfo/i observatory' observed the Low Mass X-ray Binary EXO 
0748-676^ during its Commissioning and Calibration phases for almost a half million 
seconds, spread over six satellite orbits between 2000 February 21 and 2000 April 21. 
Data were recorded with the Reflection Grating Spectrometer^(RGS) for 335,000 
seconds (data obtained with the European Photon Imaging Cameras'*'^ (EPIC) are 
available for 39,000 seconds of simultaneous exposure). During this time, a total of 28 
X-ray bursts were recijrded with the RGS, lasting a cumulative 3200 seconds. During 
the brief bursts, the neutron star outshines the accretion-generated light by an order of 
magnitude in intensit> , while the ongoing accretion ensures a continuing supply of 
heavy elements in the stellar photosphere. This makes the burst spectrum a promising 
place to detect absorption structure from a neutron star photosphere, a long standing 
goal in compact object astrophysics. Due to the long exposure, and the high efficiency 
and spectral resolving power of the RGS, this EXO 0748—676 dataset is by far the most 
sensitive to date for conducting such a search. 

Data were processed with the XMM-A^ew/on Science Analysis Software (SAS) 
that is currently available as version 5.3.3. The soft X-ray lightcurve of EXO 0748-676 
shows considerable variability*. Searching for X-ray bursts in the RGS lightcurve, we 
only considered events which conformed to the steep rise/exponential decay shape 
characteristic of type f X-ray bursts^. We identified 28 bursts in the RGS data. The 
bursts varied in peak intensity between 4 cts s ' and 12 cts s * with an average of 8.8 cts 
s’*. This represents an increase by a factor of ~15 over the quiescent levels observed in 
the periods of low acti vity. We defined the onset of a burst as the time at which the 
count rate first rises aliove the quiescent level by a factor two or more. The end, less 
well defined, occurs v'hen the count rate drops back to the local average level. The 
majority of the bursts ranged in duration from 48 to 128 seconds, with an average of 90 
seconds. Seven of the bursts were longer, with durations between 176 and 320 seconds. 
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We then extracted the first order (m = -1) RGS spectra for each burst. The 
spacecraft pointing was stable during observations on each of the separate revolutions, 
but differed between revolutions by up to 40 arcsec. To generate the average burst 
spectrum we therefore combined the data for all observations within a single revolution, 
but generated separate spectral files and response matrices for each revolution. All 
spectral fitting was pt:rformed simultaneously on these separate data sets. For ease of 
display we generated a flux spectrum of the average burst, using the effective area 
curves for each separate data set. This allowed us to combine the data from the two 
RGS instruments as v'ell. The wavelength scale is accurate to ~10 mA. The effective 
area is accurate to 5^ for all wavelengths longer than 8 A.^ Background subtraction 
was performed using the same extraction algorithms, but over the image region not 
occupied by the source. The background flux becomes a significant fraction of the total 
flux for wavelengths longer than ~32 A. We therefore only consider the wavelength 
range from 8 to 32 A in our analysis. 

To constrain the broad-band properties of the burst spectrum we examined the 
EPIC data. Pile-up during the bright bursts contaminated all but 250 seconds, or three 
bursts. The EPIC/PN spectrum of these three bursts is well fit by a blackbody, with a 
peak color temperature of ^Tbb -1.8 keV, decaying to kT^s < 1-5 keV. Since the 
spectral properties clearly evolve during the bursts, we investigated the RGS spectrum 
of the early, bright phases and the decay phases separately. We explored a variety of 
ways to subdivide the bursts, but found that the results are not sensitive to the exact 
criterion used to compile the ‘early’ and ‘late’ time burst spectrum. We therefore chose 
to divide the bursts in the simplest way, by splitting them in half by duration. The 
resulting flux spectra for the early and late phases of the averaged burst are shown in 
Figure 1. 
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As in the case cf the quiescent spectrum,® we see clear evidence for absorption 
and emission from highly ionized gas surrounding the neutron star during the bursts. O 
VII K-shell emission (consisting of « = 1-2 resonance, intercombination, and forbidden 
lines at 21.60, 21.80, and 22.10 A) is clearly detected, as is absorption by O VII 
(photoelectric absorption at 16.78 A, resonance line absorption at 21.60 A). As before, 
the fact that the O VIJ line emission is dominated by the intercombination transition 
indicates that the gas is recombining, which implies that the ionization is driven by 
photoionization, and that the electron density is relatively high (n > 10 cm' ). We see a 
clear change in the nature of the O VII spectrum as the bursts progress, with the 
emission weakening and the absorption increasing from the early to the late phases. 

The density of the plasma is sufficiently high that transient ionization effects are 
probably not important, so that a change in the relative strength of ionization and 
absorption in the sami: ion cannot be attributed to an evolving ionization balance. 
Instead, the change in the observed O VII spectrum is indicative of an overall change in 
the geometry of the source. In the quiescent spectrum, we concluded from a 
comparison of the photoionization and recombination rates, as measured by the depth of 
the edge and the luminosity of the line emission, that the O VII zone of the circumstellar 
absorber is flattened in a torus-like shape, occupying roughly 1/3 of the sky as seen 
from the neutron star." During the bursts, the weakening of the emission and the 
deepening of the absoi ption indicate a shift from a more spherical geometry at the 
beginning of the bursts toward the flattened geometry of the quiescent state at the end of 
the bursts. 

In order to develop a physically consistent model for the spectral transmission of 
the circumstellar absorber, we fit the O VII and O VIII spectra and measured the 
intensity, velocity width and Doppler-shifts of the emission lines, and the ion column 
density and velocity broadening of the absorption features for both the early and late 
phase spectra. The measurements assume an empirical continuum spectrum during the 
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bursts, constructed from a blackbody and a powerlaw, whose parameters were 
optimized by eye. O VII and O VIII absorption spectra were calculated with a spectral 
code originally developed to interpret the X-ray absorption spectra of AGN.® The model 
incorporates atomic structure and transition probabilities, and for any given ion, 
consistently accounts for the absorption in all transitions out to high principal quantum 
number and the photc'dectric absorption edge. The entire spectrum was subject to 
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absorption by a neutral medium, with equivalent hydrogen column density A^h = 1 x 10 
cm'^; the neutral absorption spectrum has strong O K absorption in the 22-24 A range, 
the shape of which we optimized to conform to the interstellar O absorption spectrum 
measured in other sources.^ 

Using the measured line intensities and optical depths, and the value of the 
photoionization parameter''^’" ^ = L\omzm%lnF^ at which the abundance of O VII peaks, 
we can parameterize ihe average properties of the circumstellar medium in terms of the 
thickness AR of the shell at radius R from the star, its average density n, and a clumping 
parameter C = </i^>/<n ". We find that the O VII zone must lie near the edge of the 
accretion disk (/? ~ 5 \ 10‘“-2 x 10*‘ cm), be relatively thin {AR ~ /?/10), highly 
clumped (C ~ 100), and fairly dense {n ~ 10'^-10‘‘* cm'^). Wavelength shifts in the O 

VII absorption feature s indicate a significant bulk outflow velocity of v ~ 5000 km s ' 
during the bursts. 

Using the measured absorption parameters for O VII and O VIII, we then 
synthesized a model for the full spectral transmission of the circumstellar absorber. We 
adopted an ionization parameter of ^ = 10 as representative, so as not to overproduce O 

VIII absorption, and added the other elements at their solar abundances, with ionization 
fractions derived from the photoionization equilibrium balance. The full set of ions 
includes the K-shell ions of C, N, O, Ne, Mg, and Si, and the L-shell ions of Fe. We 
scaled the turbulent velocity broadening of each ion with the value measured in the O 
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VII resonance line, assuming a common temperature for all ions. We adopted the 
Doppler blueshift observed in the O VII features for all ions. The resulting transmission 
model, superimposed on the optimized continuum model, is overplotted on the observed 
spectra in Figure 1. TIte apparent absence of A, 24.78 A N VII Lyoi absorption in the 
data implies a subsolar N/0 ratio in the absorbing gas. The absence in the data of 
absorption by ions that are present at higher ionization parameter — specifically the 
absence of significant Fe L absorption - implies that the circumstellar medium occupies 
only a narrow range of (fairly low) ionization parameters. 

We now examine the spectrum for any remaining structure that is not associated 
with the circumstellar absorber. In view of the noise levels, it is difficult to perform 
such a search effectively using statistical significance criteria only. We therefore appeal 
to spectroscopic consistency arguments when assessing apparent absorption features in 
the spectrum. In the i;arly time burst spectrum, the most significant modulation appears 
at 13.0 A. We also sei.’; weaker structure at 25.3, 26.3, and 26.9 A. In the late time 
spectrum, we identify significant modulations at 13.75, 25.2, and 26.4 A with weak 
features at 17.8 and 19.7 A. 

We inspected the spectra obtained by the two separate RGS spectrometers and 
find that all the features appear in both parallel spectra, but with only marginal 
significance in the case of the 26.3, 26.9 A (early times) and 17.8, 19.7 A (late times) 
features. The fact that all features appear either at early or at late times rules out the 
possibility of stationai y modulations in the spectrometer efficiency. We have also 
examined the three very deep RGS spectra of extragalactic featureless continuum 
sources (Mkn 421, 3C 273, and PKS2 155-304; Rasmussen et al., in preparation). No 
significant modulations are observed in these spectra at the positions of the features 
under discussion, nor do these spectra exhibit unexplained features at any other 
wavelength, of strengi h comparable to the ones observed here. 
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Neither do the features match the wavelengths of absorption lines expected from 
the circumstellar absorber, with the possible exception of the 26.4 A feature in the late 
phase spectrum, which is probably partly due to CVI Lyy at 26.6 A. We examined 
absorption models appropriate to higher ionization parameters (over the range ^ = 
10-100) and found that the features cannot be made consistent with circumstellar 
absorption at any ionization parameter. This explanation can be rejected on purely 
spectroscopic grounds; it requires significant velocity shifts that vary randomly between 
ions that are present at similar ionization parameters. Furthermore, fitting individual 
absorption lines to the features produced serious inconsistencies in terms of the 
predicted overall absorption structure due to any given single ion. 

We are left with the exciting possibility that the absorption features arise in the 
photosphere of the neutron star. The magnetic fields in the neutron stars in low mass X- 
ray binaries are believed to be small, so that field effects do not affect the atomic 
structure, and we can use well established atomic spectroscopy to interpret the 
wavelengths. In the only other example of nontrivial structure in the spectrum of a 
neutron star,’^ the conect spectroscopic identification, and hence the inferred redshift, 
depends critically on tlie unknown strength of the stellar magnetic field. 

We have examined all the spectra of the K-shell ions of C through Si, and find no 
multiple coincidences between observed and predicted line positions, with identical 
redshifts for all ions. However, an important clue to identifying the features lies in the 
presence of the 13.0 A line at early times and the presence of the 13.75 A line at late 
times. A solution to the Saha ionization balance for densities n ~ 10 cm' (as expected 
in a neutron star atmosphere) indicates that iron should be primarily in its H-like charge 
state at temperatures k T ~ 2 keV, while aXkT ~ 1 keV, the He-like charge state 
dominates; these temperatures roughly correspond to the observed color temperature 
early and late in the bursts, respectively. The n = 2-3 transitions in the H-like ion (the 
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analogue of the Ha Balmer line), occur at 9.518, 9.533, 9.579 and 9.675, 9.690, 9.738 
A.*^ Identifying these transitions with the feature at 13.0 A in the early phase burst 
spectrum implies a redshift of z = 0.35. The strongest n = 2—3 transitions in the He-like 
ion occur at 10.213, 10.048 A (Ehud Behar, priv. comm.; the 10.213 A transition has 
the highest oscillator strength). Identifying these transitions with the feature at 13.75 A 
in the late phase spectrum also implies a redshift of z = 0.35. The higher order series 
members will all lie at wavelengths shortward of 10 A, to which our spectrum is not 
sensitive. The Saha balance at temperatures kT> 1 keV and densities n ~ 10^^ cm'^ 
indicates that all the lighter elements should be nearly stripped, so we would not expect 
to see their K-shell absorption lines, with the possible exception of oxygen, in view of 
its relatively high abundance. Applying the same redshift to the O VIII Lya line (rest 
wavelength 18.97 A) we would expect to see a feature at 25.6 A. We speculate that the 
double XX 25.2, 26.4 A structure observed at late phases, which is centered at this 
wavelength, is in fact a self-reversed, broad O VIII Lya line, where the self-reversed 
profile is indicative of extended structure to the outer atmosphere, and possibly a slow 
outflow, such as observed in the strong UV resonance lines in massive stars with 
extended atmospheres.'"* The corresponding O VIII LyP line (21.60 A at z = 0.35) 
would be hidden in the O VII spectrum from the circumstellar medium. The low- 
significance features identified above have no obvious spectroscopic interpretation, and 
are most likely statistical fluctuations. 

We have identified three sets of redshifted transitions in iron and oxygen in the 
EXO 0748-676 spectrum, all with an implied redshift of z = 0.35. This redshift, 
combined with recent model calculations for the neutron star mass-radius relation'^ 
taken at face value, would indicate a mass of A/ - 2.0-2.2 and radius R ~ 8-10 km. 
This is consistent witli a neutron star with a birth mass near the average for non- 
accreting neutron stars (M = 1.4 Af^), that has been accreting mass at the observed rate 
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for ~ 10^ years, and is consistent with the estimated masses of other accreting neutron 

stars/*^'^ 
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Figure 1: The background-subtracted flux spectra for the early (top panel) and 
late (bottom panel) phases of the 28 bursts from EXO 0748-676 as measured 
with the RGS in the 8-32 A band. The data are plotted as the black histograms, 
with la error bars derived from counting statistics. The red line is the empirical 
continuum, with additional O VII intercombination line emission, modulated by 
absorption in photoionized circumstellar material. In red, we have labeled the 
positions of the most prominent discrete absorption lines from the circumstellar 
medium; in the He-like spectra, 'w' signifies the n = 1-2 resonance transition, 

'xy' the (unresolved) n = 1-2 intercombination transitions, while higher series 
members are marked 'KP,y', etc. Column densities in ions other than O VII have 
been normalized to the absorption measured in O VII, assuming ionization 
parameter ^=10, and solar abundances. The N VII Lya line at 24.78 A is 
overpredicted, indicating a subsolar N/0 abundance ratio. The black labels 
indicate the interstellar O 1s-2p absorption line. Blue labels indicate the 
photospheric absorption lines in Fe XXVI, XXV, and O VIII, at a redshift z = 

0.35. The data and models have been rebinned to AA, = 0.124 A, which is a 
factor of ~ 2.5 larger than the RGS instrument resolution. 
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